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Abstract 

The ionospheric plasma is described by means of the 

• • * ^ 4 * 

momentum and continuity equations for 0 , He , H and the 
energy equations for T and T. . The ion production rate profiles 

C JL 

and the photoelectron heating rates to the plasma are computed 
separately and serve as source functions in the continuity and 
energy equations respectively. The neutral atmosphere (compo- 
sition, winds, temperature) as well as the incident photon 
spectrum ate part of the inputs. Charge transfer, cojiisional 
energy loss processes, ion-neutral drag, diffusion and electron 
heat conduction are among the physical processes included in 
the calculations. Ion heat conduction is neglected. Assuming 
steady state, the equations are transformed into integral equations 
and then solved iteratively along geomagnetic field lines from 
the region of chemical equilibrium up to the equatorial plane. 

In this scheme non local heating is included in a self-consistent 
way. The lower boundary conditions are obtained by satisfying 
photochemical and local energy equilibrium. Under conditions of 
asymmetry with respect to the equatorial plane the solutions 
are carried out for both hemispheres and the upper boundary conditions 
are determined by requiring continuity of the physical parameters 
across the equatorial plane. In this way interhemispheric plasma 
transport is introduced in a natural way. For symmetric con- 
ditions the transport fluxes are assumed zero at the equator. An 
extension of this model to include time dependent phenomena is 
presented . 


iii 
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I. Introduction 


The model presented here has provided the theoretical 
framework for the interpretation of satellite and rocket data 
on the ion composition and plasma temperatures (Mayr et al. , 
1967, Mayr et al. , 1968, Brace et al . , 1969, Brinton et al . , 

1969, Mayr and Brace, 1970, Mayr et al. , 1970a, Brinton 

et al. , 1971, Mayr, 1968, Brinton et al. , 1970, Brace et al. , 

1970, Mayr et al . , 1970b) . 

In general, density and temperature gradients along 

field lines cannot be obtained from satellite data nor do 

those data provide information on the transport fluxes 

directly. For this reason it is desirable to eliminate 

uncertainties associated with the choice of upper boundary 

conditions (which require the knowledge of transport fluxes) 

as far as possible. In the interest of this objective we 

have therefore chosen to integrate the continuity, momentum, 

and energy equations along field lines up to the equatorial 

plane, thereby eliminating the boundary value problem. 

The calculations presented in this paper cover the 

F 2 regi.on and the protonosphere since satellite orbits 

rarely dip below the F 2 region. This simplifies the analysis 

in a number of respects considerably since in that case the 

model can be restricted to a three constituent plasma with 

+ + + 

the ion species, 0 , H , and H each of which could be the 
major constituent in some altitude region. Since N + , 0 2 + 
and NO + ions are observed to be always the minor constituents 
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in the topside ionosphere (Taylor et al. , 1968, Taylor, 1971b) 
they can - to a good approximation - be incorporated 
individually and a posteriori after a solution for the major 
constituents has been obtained. The distribution of these 
minor species is thereby affected by the transport velocities 
and polarization field of the major ions. Inclusion of N + 
at high latitudes where its concentration can exceed those 
of the H + and H e + ions (Taylor et aJL. , 1968) and study of 
the observed longitudinal variations in 0 ^ + (Taylor, 1971b) 
could therefore provide valuable insight into the dynamic 
properties of the ionosphere. 

In the upper ionosphere chemical and local energy 
processes do not play an important role except near the 
lower boundaries where the dominance of these processes is 
utilized to determine the lower boundary conditions. In 
the lower F and E region where transport processes can be 
neglected the continuity equations degenerate from differ- 
ential to algebraic equations, and therefore entirely dif- 
ferent solution procedures are appropriate. In this region 
the ion chemistry is evaluated separately for a multi- 
constituent plasma including 0 + , 0 2 + , N 2 + , N0 + and N + 

(Pharo et. al. , 1971, Brace et al. , 1971) . 

In describing the density and temperature structure 
of the protonosphere and its interface with the region 
non-local heating (due to escaping photoelectrons) becomes 
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signif icant . This has first been pointed out by Geisler 
and Bowhill (1965) . Since this heating mechanism depends 
on the electron density distribution in a rather complicated 
way - even with the crude assumptions that are made in this 
model - a truly simultaneous solution for the mass, momentum 
and energy equations seems impractical. For this reason the 
energy equations are solved separately for the plasma tem- 
peratures T e and T. which are then used in an iterative 
scheme with the other equations that describe the ion and 
electron densities. As will be shown, this procedure con- 
verges rapidly thus yielding fairly accurate solutions with 
only a few iterations, and in this way a physically self- 
consistent solution is obtained. 

Since tb ? model presented here treats the ionosphere- 
protonosphere system as one unit, including the effects of 
charge exchange and non-local heating, without any ad hoc 
upper boundary condition, the effects of the ionosphere- 
protonosphere coupling are obtained in a natural way. The 
significance of the charge exchange reaction for this 
coupling was first pointed out by Hanson and Ortenburger (1961) 

Finally, the model is designed to describe primarily 
steady state conditions. However, it has been shown, that 
this limitation can be removed. A solution for the diurnal 
variations in the ion composition, involving plasma transport 
between F n region and protonosphere, has been presented by 
Mayr (1968) and is briefly discussed in Section IV C. 
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I 1 • Definition and Organization of the Model 

In Table I the physical processes, considered in the 
continuity, momentum, and energy equations, are summarized 
since they have been discussed in several publications. 

Figure 1 shows the organization of the model in the 
form of a block diagram. Blocks and circles with heavy 
lines represent processes and internally computed " inputs'* 
respectively, while circles with thin lines represent the 
true inputs. 

Adopting a photon distribution based on Hinteregger 
et al. (1965) and a model for the gas temperature and neutral 
composition (Jacchia, 1964) , photon absorption is considered 
and the ionization rates for 0 + and He + are computed. At 
lower heights where chemical reactions dominate, the con- 
centrations of 0 + and He + are determined from photochemical 
equilibrium. This is the case for altitudes up to 180 km 
for 0 + and up to 400 km for He + . For H + the charge exchange 
equilibrium relation involving the densities of 0 + , O and 
H is employed up to 400 km. These chemical equilibrium 
densities se/ve as the lower boundary conditions for the 

continuity equations of those constituents. 

t 4* *4* a 

The charge transfer reaction O + H -► H + O is the 
primary source for protons and accordingly the proton content 
of the magnetosphere is significantly influenced by the con- 
centration of atomic hydrogen at exospheric levels. It had 
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been postulated for several reasons that bear upon the 
energetics (Brace et al. , 1967) and ion composition (Mayr 
et ai . , 1967, Brinton et al., 1969), that the hydrogen con- 

5 

centration should be of the order of 1 x 10 /cc at exospheric 
heights during minimum solar activity. This value which has 
been used in our model is significantly higher than those 
in the atmospheric models of Jacchia (1964) and CIRA (1965) , 
but it is consistent with recent results (Brinton and Mayr, 
1971a, b) in which for the first time the concentration of 
atomic hydrogen was deduced from in situ measurements of ion 
and neutral constituents at around 350 km. 

A meridional wind field (described in terms of spherical 
harmonics by Volland and Mayr, 1963) is adopted as input for 
the calculation of the neutral-ion drag interaction parallel 
to the magnetic field. It is assumed that the wind field 
is defined in the reference frame of geographic latitudes 
which, to a first order approximation, neglects the feedback 
effect on the winds due to ion-neutral drag, which is strongly 
controlled by the geomagnetic field. 

The wind field and the plasma distribution are symmetrical 
with respect to the magnetic equator only if the latitudinal 
pressure bulge in the neutral atmosphere occurs there. Linder 
this particular condition the transport velocities are zero 
at the equator. This condition determines the solution of 
the continuity equations uniquely. When the bulge is not at 
the geomagnetic equator, the wind velocities are different 
in the northern and southern hemispheres along a particular 
field line, thus inducing asymmetries in the plasma distribution 
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with respect to the magnetic equator. Such asymmetries 
are associated with plasma transport across the equator 
which has to be considered in a self-consistent way in 
order to find unique solutions. 

The scheme in which these transport fluxes are incor- 
porated is briefly the following: For the hemisphere in 

which the plasma densities are higher on outward flow of 

+ 

the major ion component (normally H at the equator) is 
introduced and its equatorial density is computed for a 
number of flux values. For the same value inward flows 
into the field tube of the other hemisphere are then used 
to compute the corresponding equatorial densities. From 
the condition that equatorial densities must be equal, the 
inter-hemispheric transport flux for the major constituent 
is then determined. Similar procedures are applied for 
the minor constituents so that eventually the various 
transport fluxes satisfy the upper boundary condition of 
continuous density transition across the equator for all 
ion species. 

The plasma temperatures T g and T^ are initially assumed 
equal to the gas temperature to start the solution for the 
equations of mass and momentum conservation. The resulting 
electron density distribution, N 0 , is then used to compute 
the non-local heating rate due to escaping photoelectrons. 
The calculation of the local and non-local heating rates 
has been based on the method of Fontheim et al. (1960) . 
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With these heating rates and the electron-ion, electron- 
neutral, and ion-neutral loss rates, the energy equations are 
solved for the electron and ion temperatures. The boundary 
conditions are obtained by assuming local energy equilibrium 
(neglect of heat conduction) at lower altitudes and zero heat 
flux at the equator under conditions where the plasma dis- 
tribution can be assumed to be symmetrical with respect to 
the magnetic equator. Under conditions of plasma density 
asymmetry a symmetric temperature distribution has been 
adopted to study the effects of inter-hemispheric mass 
transport. The effects of inter-hemispheric energy transport 
have not been treated so far but they can, of course, be in- 
corporated into the model. 

The resulting temperature distributions are then fed 
back into the momentum and continuity equations, and the 
iteration (including the electron density dependent heating 
rates) is continued until the relative differences between 
successive density and temperature values become smaller 
than a few percent throughout the entire altitude range. 

This procedure converges very rapidly, and it is found that 
four to five iterations are usually sufficient to produce 
fairly accurate solutions. 

In this iterative cycle the calculation of the heating 
races must be performed three to four times. The computation 
time for this calculation exceeds by far the total computation 
time for the solution of the continuity and energy equations. 
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In the latitudinal temperature study reported in Section 
IV B a simplified heat input model has therefore been used. 

By assuming a monoenergetic beam of photoelectrons of 
typically 10 ev (corresponding to the maximum in the photo- 
electron spectrum (Fontheim et al. , 1968) ) and with zero pitch 
angle the non-local heat input is calculated from the Coulomb 
loss rate given by Butler and Buckingham (1962) and Dalgarno 
et al . (1963) . The local heat input is assumed to be pro- 

portional to the atmospheric density with the proportionality 
factor as an input parameter based on the heating rates of 
Fontheim et al. (1968) . 

Although thermal diffusion (Walker, 1967) has not been 
considered in the model calculations presented in Section IV, 
it is included in the present version of the model. 

III. Solution Proi ~‘.dure 

The procedure for solving the momentum and continuity 
equations is outlined for a two constituent plasma of 0 + 
and H + ions. The continuity and momentum equations have 
the form 


Q - L [0 + ] - B jL 

a S 



3 [ 0 + ] 


(i) 
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3 [H + ] 
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where N = [0 + ] + [H + ] 

C- 

tO + ] , [H + ] = concentrations of 0 + , H + 

s = distance along a field line 
Q = photoionization rate for 0 + 

L = o + loss rate for charge transfer and ion molecule 
reactions with 0 2 and N 2 (function of the concentra 
tions of H, 0^/ and N 2 ) 

B = it (1 + 3 COS 2 e) 1/2 
r 

M = magnetic dipole moment of the earth 
r = geocentric radius 
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0 = geomagnetic (dipole) latitude 
t = universal time 

+ + 

V 0 +, V H+ = transport velocities of 0 , H 

K l' K 2 = c ^ ar 9 e transfer coefficients 

v = meridional wind component parallel to the magnetic 
field 

[N] = concentration of the neutral atmosphere 

Q | 

OH = drag coefficient for momentum transfer between 
0 + and H + 

0^+ XT/ “ drag coefficient for momentum transfer between 0 + , 

ON H N 3 

• H and N 

T , = electron and ion temperatures 

m^t, m u + = masses of 0 H , H + ions 

U n. 

g^ = gravitational acceleration parallel to the magnetic 
field 

k - Boltzmann's constant 

After formally integrating equations (1) and (2) one 
obtains the following expressions for the transport velocities 


V 
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B 
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where j Q + and j H + refer to particle fluxes and the subscript 
e denotes the values of the functions at the equator. 

If the left hand sides of Equations (3) and (4) are 
formally denoted by [0 + ] D Q + ([0 + ], [H + ] , j Q +, j H +) an< ^ 

[H ] D h + ( [0 ] , [H ] , j Q +, j H +) / the momentum equations take 
the form of homogeneous equations 


k T. 
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B s 


3T B T 


B s 
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( 8 ) 


which can be integrated to describe the density distributions 
in the form 


[0 + ] = [ O + ] 0 exp 
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where [0 + ] Q and [H + ] q are the lower boundary conditions. Since 
Dq+ and D h + are integral functions of the densities [0 + ] 
and [H + ] , the equations (9) and (10) are integral equations. 

The solution of these equations is accomplished by means 
of an iterative procedure in which the functions D Q + and 
D„+ are treated as perturbations in the following scheme. 

For the kth iteration the perturbation functions are 


D o + 


e k D o^ 


d h + 


,-k-l 

e k °H + 


where is a parameter that is assumed to vary between 


0 i e k-i < e k < e k+i i 


For the initial solution we assume = 0 which is equivalent 
to diffusive equilibrium. In subsequent iterations is 
gradually increased to allow the perturbation function to 
become effective slowly. This procedure assures the applica- 
bility of the perturbation method which demands small effects 
from the "perturbing" function. A successive application of 
this procedure leads to a convergence to the solution of 
equations (9) and (10) when approaches 1. 

These solution procedures can be directly extended to 
incorporate He + ions. Figure 2 shows the results of a sample 
calculation of composition profiles including He + * 
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For the time-dependent case (-^ t 0) this scheme is 
extended by substituting for the kth iteration 


a [o + r 

~TT~ ) 


e k 


a [0 + f 

r) t , 


k-1 




W-l 


In the initial solution it is assumed that = 0 corresponding 

to steady state. The density distributions are computed in 

time intervals of 1 hour to generate the quasi-stationary 

diurnal variation under the influence of the diurnal variations 
+ + 

of the O and H boundary conditions, the solar flux, and the 
neutral atmosphere. In subsequent iterations is then in- 
creased toward 1 to produce finally the time-dependent solution. 
The energy equations for T and T. have the form 


| K e B -1 (I 

7 e 3 s \B 


a 

a s 


T? e 2 -Qe -K ei n e -T,) -K en (T e -T n ), 


( 12 ) 


Kei (T e -Tj) =K ln (Ti -T n ), 


(13) 


where 


= heat conduction coefficient for electrons 

Q e = heating rate for electrons (dependent on electron 
and neutral densities) 
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K ei/ K en' K in = rate coefficients for energy loss between 

the various electron, ion and neutral con- 
stituents thus depending on their densities. 


Equation (12) can be formally integrated to yield 


T, <s, =■ {IV * ^ 




B • 

K K ' ' n 


r s 




B (Qp 


(14) 


,2 7 


“ K „ (T p -T n ))ds 


d s' 


with 

T = lower boundary condition, computed from local energy 
equilibrium 

= heat flux across the equator normally assumed to be zero. 

The integral equation (14) is solved iteratively by 
increasing e from 0 to 1 , in a similar way as has been des- 
cribed for the solution of the ion continuity equations. 

In the following a few model studies are presented 
which demonstrate some of the dynamic and energetic charac- 
teristics of the upper ionosphere. 

IV. Discussion 

A) The Effects of Interhemispheric Plasma Transport 
Bearing on Longitudinal Variations in the Ionosphere. 

For illustration we shall consider wind effects originally 
introduced by King and Kohl (1965) during equinox to produce 
a condition of asymmetry. It is assumed that the meridional 
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wind velocity is zero at the subsolar point which in turn 
occurs at the geographic equator. Following Volland and 
Mayr (1968) the wind field is then described in terms of 
spherical harmonics of low degrees 


V-V P COS " + V 

10 10 21 


* 21 

cos 


n 


CO s 


( T - 


(15) 


where 

P 1A = sin 6 
10 g 

P A , = 4 sin 2 6 
212 g 

0 = geographic latitude 

9 

t = t + A/w (local time) 
t = universal time 

u) = 2 tt/ t q (angular frequency with t q the period of one day 
A = geographic longitude 

t 2 = 1400 local time (time of the observed diurnal 
maximum in the atmospheric density) 

v = meridional wind velocity 

Assuming that the altitudinal variations of v within the 
F 2 region are small compared with the variations in the neutral 
density and therefore negligible. Equation (15) is defined with 
two parameters v^q and ^21' according to Volland and Mayr 

(1968) are assumed to be of the order of v^q = 10 m/sec and 


50 m/sec. 
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It is the wind component parallel to the magnetic field 

which enters into the ion continuity equations. The wind 

field in (15) must therefore be described as a function of 

geomagnetic latitude. This is straightforward for the two 

longitudes in which the earth's axis of rotation, the magnetic 

dipole vector, and the earth-sun-line lie in the same plane. 

At American longitudes the magnetic equator is south of the 

geographic equator, and thus the wind field can be described 

by substituting in Equation (15) 0=9+ 12°, where 0 is 

y 

the geomagnetic latitude. This then exhibits explicitly the 
asymmetry of the wind field with respect to the magnetic 
equator. The downward wind component in the southern hemisphere 
is larger by almost a factor of two than in the northern hemi- 
sphere at about 30° geomagnetic latitude. 

With this wind field the ion composition has been computed 
along the field 1> with L = 1.4 which corresponds to about 
+30° latitude at V xegion heights. Figure 3 shows the 

^ A. 

results, on the left hand side for the northern hemisphere 
and on the right hand side for the southern hemisphere. It 
is apparent that in the northern hemisphere (left) the F 2 
peak density is larger and occurs at a higher altitude than 
in the southern hemisphere (right) , thus reflecting the aif- 

i 

ference in the wind velocities. The larger downward wind 
component thereby drags the F 2 region ionization (0 + ) to 
lower heights where the increase in the recombination rate 
causes the ion density to decrease. 
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The dotted lines show the ion distributions for 0+ and 
H which were computed assuming zero plasma transport across 
the equator. As a result it is seen in Figure 3 that the 
equatorial H + concentration is by about a factor of five 
larger in the northern hemisphere than in the southern hemi- 
sphere. Such a discontinuity across the equatorial plane 
is of course unrealistic. Therefore H + fluxes are introduced 
across the equatorial plane from the northern to the southern 

hemisphere in a scheme that has been described in Section II. 

. + 7 2 

In this case an H flux of 5.4 x 10 /cm sec is required 

to match the H + densities from both hemispheres at the equator. 
The resulting electron and H + distributions are shown in solid 
lines. The 0 + distribution is represented by the dashed line 
which coincides with the '.olid line at lower altitudes. It 
is apparent from this that the interhemispheric plasma trans- 
port, induced by the asymmetry in the thermospheric wind 
field, decreases the protonospheric density by about a factor 
of two in the northern hemisphere thus raising the 0 h -H + 
transition level by 90 km while it increases the proton 

density in the southern hemisphere and accordingly lowers 
-f + 

the O -H transition level by 70 km. The effect of this inter- 
hemispheric plasma transport is also seen in the scale height 
of the 0 + concentration within the F 2 region. At these alti- 

_l_ 8 2 

tudes the H flux is about 10 /cm sec an increase by a factor 
of 2 from the equatorial value due to the convergence of the 
field tube. The H + -0 + drag resulting from this upward proton 
flux in the northern hemisphere causes an increase in the 0 -1 ’ 
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. -f- 

scale height above the region where H is in chemical equi- 
librium ( up to about 500 km) . Conversely in the southern 

hemisphere the downward flux decreases the 0 + scale height. 

+ 

Below 500 km where H becomes chemically controlled, the 

ionization transport is taken over from H + to 0 + ions thus 
+ 8 

producing O fluxes of the order of 10 above the F 2 peak. 

These fluxes therefore eventually influence the height and 
density of the F 2 peak, though only slightly. 

Up to this point only interhemispheric plasma transport 
by protons has been considered. Under this condition the 0 + 
density is still discontinuous across the equator as is apparent 
from Figure 3 in which the 0 + density distribution is shown 
as a dashed line above 1400 km. It is apparent that, with 
the H + flux alone, the 0 + densities from both hemispheres differ 

by more than a factor of 10 at the equator. However, after 

. + + 
the introduction of 0 fluxes, the density continuity in O 

5 2 

is established with a flux of 5.1 x 10 /cm sec. The resulting 
density distributions are shown in solid lines in Figure 3. 

The effect of this 0 + flux is negligible for the H + distribution, 
and it affects 0 + only above 1600 km as is evident from the 
deviation between the dashed and the solid lines in Figure 3. 

It is apparent from this that in the northern hemisphere, 
where the F 2 region densities are higher, the 0 + scale height 
decreases significantly above 2000 km, while in the southern 
hemisphere, where 0 + is supplied, the concentration of 0 + 
even increases with height. This somewhat peculiar distribution 
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thereby indicates that the minimum in the 0 + density is 
shifted from the equator toward the southern hemisphere 
which is quite plausible considering that the F 2 region 
ionization is also lower there. 

At other longitudes, and in particular at European and 
Asian longitudes, the asymmetry in the ion density distribu- 
tions between the southern and northern hemispheres would be 
opposite to the case discussed here. As has been shown by 
Brinton et al. , (1970) , this effect can account for signifi- 

cant longitudinal variations in the 0 + -H + transition level. 
Taylor et al. , (1970) and Taylor (1971a) have presented 

comprehensive evidence for the solar geomagnetic control of 
the ion composition from OGO 4 and 6 ion mass spectrometer 
measurements thus demonstrating that the longitudinal vari- 
ations - which bear upon the effects discussed here - are 
very significant and are comparable in magnitude with the 
latitudinal and local time variations in the ion composition 
(Taylor, 1971b) . 

B) Latitudinal Electron Temperature Variation 

Equations (5) , (6) , (9) , (10) , (12) and (13) 

have been solved for several latitudes under equinox condi- 
tions (zero flux across the equatorial plane) to study the 
latitudinal variations in the electron temperature (see 
for detail.? also Mayr et al . , 1968) . In Figure 4 the electron 
temperature is shown as a function of geomagnetic latitude 
for several altitudes. For comparison, the observed lati- 
tudinal temperature distribution at 1000 km (Bi e eit al . , 
1967) is also presented. The agreement between the measured 
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and calculated temperatures is excellent up to 60° geo- 
magnetic latitude. In particular, the calculated profiles 
clearly show the equatorial trough and the maximum at mid- 
latitude. Above 60°, at the plasmapause and beyond, H + is 
depleted by electric field convection and polar wind escape. 

This mechanism is not included in our equations, and there- 
fore the disagreement with experiment above 60° is to be ex- 
pected. 

To provide some insight into the physical mechanisms 
which determine the downward heat flow at 1000 km the various 
integrated energy rates are shown in Figure 5. The dashed 
line shows the energy input due to fast electrons, the dash- 
dotted line shows the energy loss between plasma and neutral 
atmosphere, and the solid line represents the difference be- 
tween input and loss, which represents the heat flux at 1000 km. 

The latitudinal variation of the energy input above 
1000 km depends strongly on the rate at which the photoelec- 
trons lose their energy at lower altitudes. Near the equator, 
the electron density is high, and almost all of the energy is 
absorbed below 1000 km. Thus only a small amount of energy 
can heat the ionosphere above this height. With increasing 
latitude the electron density decreases r and an increasingly 

t 

larger fraction of the energy becomes available as input at 
higher altitudes. Above 40° the latitudinal variation in 
the electron density is small, and therefore the energy input 
above 1000 km stays constant. 
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The total rate of energy loss due to collisions in 
field tubes above 1000 km (also shown in Figure 5) depends 
proportionally on the ion content and on the density of the 
neutral atmosphere (primarily neutral hydrogen) • Addition- 
ally, the loss rate is a function of the ion and neutral 
species involved in this loss process. Hydrogen ions lose 
their energy most effectively to hydrogen atoms, and similarly 
oxygen ions cool most rapidly to neutral oxygen but at a 
rate 4 times lower than the loss rate between H + and H (Brace 
et al . , 1965). At the equator, the loss rate is small because 
of the small field tubes and their small particle population. 
With increasing latitude the field tube volume increases and 
so does the total ion content. The result is an increase of 
the energy loss rate. Between 20° and 40° latitude the ion 
content increases further, and there is a small increase in 
the population of neutral particles. But simultaneously 
the mean ion mass increases thus decreasing the effectiveness 
of the cooling process. The net result is that the loss 
rate changes little with latitude. At latitudes above 40° 
the field tube volumes strongly increase and therefore the 
contact between plasma ai i neutral atmosphere involves 
significantly larger particle populations. The net effect 
is that the integrated loss rate rapidly increases at higher 
latitudes . 

The variation of the heat flux (solid line in Figure 5) 
reflects the very different behavior of input and loss rates 
with latitudes. Near the equator both energy input and loss 
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rates are low. The loss exceeds the gain, and this requires 
a negative temperature gradient along field lines that con- 
tributes to the low equatorial temperature. Up to middle 
latitudes the energy input increases while the loss function 
changes little. Therefore the heat flux increases as does 
the temperature. Above 40° the energy loss starts to increase 
and the input levels off. Thus the heat flux forms a maximum 
at 50° where T g is also highest. Beyond this latitude the 
energy loss rate rises rapidly with latitude. This leads to 
a rapid decrease of the heat flux and to a similar decrease 
of the temperature at high latitudes. 

Figure 4 shows that at high latitudes (above 60°) , the 
calculated temperatures decrease significantly below the 
measured ones. The reason for this discrepancy may be that 
some of the assumptions in our model are not valid at high 
latitudes. For example, electric fields which would give 
rise to Joule heating may play an important role outside the 
plasmasphere; yet their effects have not been included nor 
have we considered the effects of precipitating electrons 
in the present model (Paulikas, 1971;, Fontheim et al., 1971a; 
Fontheim jit ad., 1971b) . 

The ion composition which was simultaneously computed 
with the electron temperature was used to determine the mean 
ion mass, m + , at 1000 km. It is shown in Figure 6 in com- 
parison with results deduced from OGO 2 (Taylor et al. , 1968) . 
The agreement is good up to middle latitudes. However above 
50°, where the discrepancy in the temperature becomes apparent, 
the observed mean ion mass increases to its maximum value of 
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16 AMU, whereas the computed values of m + decrease to 3 AMU. 
This means that in our model the electron density scale height 
is too high by a factor of 5 and consequently the total ion 
content above 1000 km is also too high. This discrepancy 
has been interpreted as due to a continuous upward flux of 
protons that increases the mean ion mass and raises the 
0 + -H + transition level (Mayr et al. , 1967) thus masking 
at these high latitudes the diurnal redistribution between 
ionosphere and protonosphere (Mayr, 1968) . These fluxes 
which have also been inferred from the ion distributions 
across the plasmapause (Taylor et al. , 1969, Mayr et al . , 
1970a, Brinton et al. , 1971) are associated with electric 
field convection (Nishida, 1966), and polar wind escape 
(Banks and Holzer, 1968) . 

As shown earlier, the integrated collisional energy 
loss rate is very sensitive to the total plasma content in 
the field tube and to the mean ion mass. A reduction of the 
proton content would lead to a reduction in the integrated 
energy loss rate, and this, in turn, would tend to increase 
the temperature. A significant reduction in the proton 
content would also cause the protonosphere to become trans- 
parent to photoelectrons and as a consequence the heat input 
into the p_otonosphere would decrease as well. Most of the 
photoelectron energy could therefore be dissipated in the 
conjugate hemisphere thus heating primarily the plasma below 
1000 km which in turn would reflect upon the temperature at 
this height level. 
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Hence, the discrepancies between the calculated and 
measured composition and temperatures at high latitudes could - 
though in a rather complicated manner - be related to the same 
physical mechanism, the loss of plasma outside the plasma- 
sphere . 

C) Discussion of Time-Dependent Solutions 

The time-dependent equations are solved iteratively 
as has been discussed in Section III (for details see Mayr, 
1968) . In this way the problem is again reduced to a boundary 
value problem. At the lower boundary (300 km) the temporal 
variation of [O ] is assumed to be of the form 


[0 + 3 = [O + ] 0 + [0 + ], sinlUt, (16) 

T o 

where t q is the period of one day. H + is assumed to be in 
chemical equilibrium through charge exchange, and its density 
at the boundary is therefore given by the expression 


[H *] = * Lqp-M 
8 [ 0 ] 


(17) 


The neutral densities [0] and [H] at 300 km are assumed to 
be given by an expression of the form (16) and so are the 
neutral and charged particle temperatures. The ion and 
electron temperatures are assumed to be equal and to be 
constant with altitude. 

Calculations have been performed along a midlatitude 
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Table II lists the values of the input parameters. In 
Figures 7 through 9 the time-dependent solutions (solid 
curves) for 0 + , H + , and electrons a.’-e compared with the 
quasi-stationary solutions O/Bt = 0) (dashed curves). It 
is immediately evident that the differences between the time- 
dependent and stationary solutions are significant (especially 
at high altitudes) , and thus the steady state assumption is 
not always appropriate in describing the topside ionosphere. 

The 0 concentration (Figure 7) exhibits a noon maximum 
at all altitudes. Since the same lower boundary condition 
is being used for the steady state and time -dependent cases, 
the solutions hardly differ at lower altitudes. With increas- 
ing altitude the effect of the boundary condition becomes 
less important, and therefore the time-dependent and stationary 
solutions start to deviate. 

The H + density (Figure 8) shows a significantly more 
complex behavior. At 300 km its peak is in phase with 0 + , 
and its minimum stretches over several hours during the 
night. This is a consequence of the ion chemistry that in- 
volves variations of 0 + , O, and H. At 500 km, however, the 
temporal behavior of H + shows a radically different picture. 

A broad minimum at daytime is followed by a sharp maximum 
during the night. The reason for this change can again be 
found in the ion chemistry which at this altitude involves 
significantly different [0]/[H] ratios. In addition it is 

evident from the difference between the stationary and time 

. H- 

dependent solutions that downward fluxes of H enhance the 
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proton concentration at night by almost a factor of two. 

At 700 km the temporal variations of H result from the 
interaction of several mechanisms. At this altitude the 
transition from the region to the protonosphere occurs 
which means for H + that transport processes start to dominate 
chemical processes. This is apparent from the large dif- 
ferences between the stationary and time dependent solutions 
which show that during daytime upward fluxes decrease the 
H + concentration while at night downward fluxes enhance the 
H + concentration. These features prevail up to about 2500 km 

* “f- 

where the H density is decreased during the day and enhanced 
at night such that the amplitude of the temporal variation 
is lowered by a factor of two and the phase is shifted by 

-f + 

almost three hours. The diurnal behavior of O and H dis- 
cussed above is in essential agreement with the measurements 
of Brinton et al., (1969). 

The dynamic effects are exhibited best by the electron 
density distribution (Figure 9) . It shows that with increas- 
ing altitude the amplitude decreases while density minima 
and maxima shift in time. These effects occur at night above 
500 km and at daytime above 700 km, the reason being that 
the 0 + -H + transition region is at a lower altitude at night 
than it is during daytime. 

In Figure 10 the altitude of the 0 + -H + transition level 
in the time-dependent case is compared with that for the 
steady state case. The difference is very pronounced show- 
ing that dynamic effects raise the transition level during 
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the day by 150 km and lower it at night by 100 km. Thus 
the amplitude of the variation of the transition level 
increases as the result of dynamic effects, which give 
rise to upward fluxes of plasma during the day and down- 
ward fluxes at night. The nighttime downward flux has also 

been predicted by Nagy et a_l. , (1968). In general, the 

. + 

daytime upward flux of H above the transition level is 
sustained by an upward flux of 0 + between the F ^ region and 
the transition level. At night, of course, the flows are 
reversed. This diurnal variation of the o"’’ flux has been 
observed b^ Evans (1971) and by Carpenter and Bowhill (1971). 

V . Summary and Conclusion 

The icnosphere-protonosphere system is governed by 
many interlocking mechanisms. Some of these are chemical 
reactions (including ionization and recombination), trans- 
port of both mass and energy, energy input, and energy loss. 

In addition, the neutral atmosphere interacts with the ion- 
osphere and is in turn modified by it. The various mechanisms 
are determined by cross sections, rate, and transport coef- 
ficients most of which are not very accurately known. In 
setting up the equations describing the ionosphere-proton- 
osphere system it is, therefore, neither practical nor 
physically meaningful to write them in their most general 
form. Fortunately, not all of the mechanisms are of equal 
importance everywhere. Significant simplifications can, 
therefore, be achieved by restricting a given model to a 
certain range in altitude or latitude or both. 
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Tne equations discussed above apply to the ionosphere 
at low and midlatitudes (inside the plasmasphere) above the 

F region. Therefore the equations involve only three con- 

+ + + 
stituents, 0 , H and He 

The emphasis in this treatment is on the dynamic effects. 
It is, therefore, of special importance to eliminate, as far 
as possible, any arbitrariness, arising out of the uncer- 
tainties in the boundary conditions. The lower boundary 
conditions are determined by utilizing the predominance of 
local chemical processes below the F region. The equations 
aro integrated from both hemispheres up to the magnetic 
equatorial plane. If the neutral wind fields are asymmetric 
with respect to the equatorial plane, they give rise to 
asymmetric ion density profiles. In that case interhemispheric 
ion fluxes are introduced in a self-consistent way by re- 
quiring continuity of ail physical parameters across the 
equatorial plane. These continuity requirements are the 
upper boundary conditions . 

The energy equations for T g and T^ are solved separately 
in an iteration in series with the continuity and momentum 
equations. This scheme facilitates the inclusion of n^;.- 
local heating which depends on the density distribution of 
the ambient plasma. 

Some of the conclusions drawn from the case studies 
presented here are summarized below: 

1) Even under steady state conditions the ion transport 
fluxes induced by wind generated asymmetries in 
the plasma population are significant and thus affect 
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considerably the o + -H* transition level and the 
scale heights of all ion species. In general, 
the assumption of diffusive equilibrium does there- 
fore not hold, not even for the major ion consti- 
tuents . 

2) The latitudinal variation of electron temperatures 
at 1000 km is shown to be related to the processes 
of non-local heating and collisional energy loss 
both of which depend in a complex way on the electron 
and ion distributions within each particular field 
tube . 

3) If the time dependence is taken into account, the 

diurnal plasma exchange between the F 2 region and 

the protonosphere is shown to be significant thus 

affecting considerably the diurnal amplitude of 
+ + 

both the 0 -H transition level and the protonosphenc 
density content. 
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Table I 

Physical Processes 


1) Photoionization 
0 + hu ■+ 0 + + e 
He + ho ■* He + + e 

2) Charge transfer and ion molecule reactions 
0 + + 0 - -*• 0 + + 0 

0 + + N_ -*■ N0 + + N 

H + + 0 + 0 + + H 

He + +*N ■* N + + He 

3) Drag (diffusive) interactions between the various ion species 

4) Drag interaction between neutral and ion species, considering 
the meridional wind component along field lines 

5) Collisional energy loss processes between the various 
charged and neutral constituents 

6) Thermal heat conduction in the electron component; ion heat 
conduction is neglected 

7) Local and non-local heating 

8) Interhemispheric transport of plasma under conditions of 
seasonal asymmetry 
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Tabl 

List of Input Parameters 

[0 + ] q = 2.0 x lO^cm ^ 

[0] q = 1.5 x lO^cm 2 

[H] = 2.0 x 10 5 cm" 3 

T. = 2.0 x 10 3 °K 

T n0 = 1.0 x 10 3 °K 


II 


for Time-Dependent Case 

[0 + ] i = 5.0 x lO^cm 2 
[0] ^ = 1 . 0 x lO^cm 2 

[H] ^ = 2 . 0 x lO^cm ^ 

T 1 = 5.0 x 10 2 °K 

T nl = 2,0 x 1o2 ° k 
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Figure 1. Organization of the ionosphere model. 
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Figure 2. Calculated composition profiles of 0 + , H + , He + , and electrons. 
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GEOMAGNETIC LATITUDE 

Figure 4. Computed latitudinal electron temperature distribution for various altitudes 
(solid curves). For comparison the electron temperature at 1000 km from 
Explorer XXII, vernal equinox 1965, (Brace et al., 1967) is shown in the 
dashed curve. 
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Energy input above 1000 km 

Energy Loss above 1000 km 

Heat Conduction flux at 1000 km 
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GEOMAGNETIC LATITUDE 

Figure 5. Energy rates integrated over the field tube volumes above 1000 km (with base 
area of 1 cm^) as a function of the latitude a: the base of the field line. 
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Figure 6. Computed mean ion mass at 1000 km. For comparison the mean ion mass 
deduced for 1000 km from OGO 2 results (Taylor et al., 1968) is shown for 
dusk, October 1965. 
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Figure 9. Computed temporal electron density distribution at various altitudes. The 
solid curves represent the time-dependent solution and the dashed curves 
the quasi-stationary solution (d / Jt = 0). 





